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amplitude  of  initial  disturbance 
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skin  friction  coefficient,  2t  /pU 

w 
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aerodynamic  lift  coefficient,  L/^pU  x  Area 

compressor  or  turbine  blade  chord  length 

shape  factor,  6  /0 

blade  incidence  angle 

surface  roughness  height 

equivalent  sand  grain  roughness  size 

length  of  laminar  region  in  separation  bubble  up  to  transition, 

x  —  x 
t  s 

length  of  transition  region  in  separation  bubble  up  to  the 
turbulent  reattachment  point ,  x^,  -  xt 
lift  force 

pressure  gradient  parameter,  (0  /v)  dU/dx 
upstream  stagnation  pressure 
chord  Reynolds  number,  Uqc/v 


crit  instability  Reynolds  number,  Ux^/v 

boundary  layer  thickness  Reynolds  number,  U6/v 

* 

displacement  thickness  Reynolds  number,  Uo  /v 
momentum  thickness  Reynolds  number,  U0/v 
Tu  freestream  turbulence  intensity  parameter. 


U 

local  x-direction  boundary  layer  velocity 
local  x-direction  velocity  fluctuation,  rms 
boundary  layer  velocity  at  top  of  roughness  element 
local  freestream  velocity 


incoming  freestream  velocity 

\ 

freestream  velocity  at  the  separation  point 
general  streamwise  location 

streamwise  location  of  instability  point  of  neutral  disturbances 

streamwise  location  of  incipient  laminar  separation 

streamwise  location  of  transition 

streamwise  location  of  fully  turbulent  flow 

Greek  Symbols 

boundary  layer  thickness 
displacement  thickness 

intermit tency  factor  for  developing  turbulence 

2 

Pohlhausen  pressure  gradient  parameter,  (6  /v)  dU/dx 
momentum  boundary  layer  thickness 

momentum  boundary  layer  thickness  at  the  incipient  laminar 

separation  point 

wall  shear  stress,  y(du/dy)^_Q 

kinematic  fluid  viscosity 

fluid  density 

absolute  viscosity 


1.  INTRODUCTION 


This  paper  summarizes  a  review  of  the  literature  on  laminar-to- 
turbulent  transition  phenomena  in  incompressible  flow.  Emphasis  has  been 
placed  on  the  transition  process  in  subsonic  axial-flow  turbomachine 
blade  rows,  particularly  compressor  blade  rows. 

Accurately  predicting  transition  in  turbomachine  blade  rows  is 
important.  Increasing  emphasis  has  been  placed  on  modeling  the  loss 
characteristics  of  axial-flow  compressors;  these  losses  are  correlated 
to  the  momentum  thickness  of  the  boundary  layer  at  the  trailing  edge  of 
the  blades.  Since  the  trailing  edge  boundary  layer  is  often  turbulent, 
the  momentum  thickness  depends  on  location  of  the  laminar-to-turbulent 
transition  region.  Also,  the  performance  of  compressors  is  adversely 
affected  by  the  presence  of  laminar  separation  on  the  airfoil  suction 
surface.  Breakdown  of  the  laminar  boundary  layer  to  turbulence  prevents 
laminar  separation  and  the  accompanying  performance  deterioration.  Pre¬ 
diction  of  the  transition  location,  therefore,  is  an  important  requirement 
for  compressor  performance  modeling. 

Although  the  transition  process  has  been  studied  for  many  years,  the 
exact  mechanism  of  transition  is  not  fully  understood.  In  fact,  no  suc¬ 
cessful  analytical  theory  of  boundary  layer  transition  has  been  developed. 
Transition  continues  to  elude  theoretical  prediction  in  terms  of  its  ini¬ 
tiation,  duration,  and  forcing  mechanism.  On  the  macroscopic  level,  a 
number  of  different  criteria  are  used  to  determine  the  onset  of  transition 
experimentally. 

Two  principal  approaches  are  used  to  provide  analytical  transition 
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prediction.  The  first  approach  is  linear  stability  theory,  which  predicts 
transition  through  the  growth  of  unstable  infinitesimal  disturbances  in 
the  boundary  layer.  The  second  approach  is  turbulence  modeling,  which 
analyzes  the  developing  presence  of  unstable  waves  through  turbulent 
boundary  layer  equations. 

This  paper  first  reviews  current  analytical  studies  regarding  the 
transition  process,  including  relevant  analysis  of  laminar  separation, 
and  then  covers  analytical  and  experimental  studies  concerning  parametric 
effects  on  transition  pertinent  to  axial-flow  turbomachine  blade  rows. 
These  include  the  effects  of  pressure  gradient,  freestream  turbulence, 
surface  roughness,  heat  transfer,  and  combined  effects.  Finally,  the  work 
of  various  authors  who  have  established  data  correlations  for  transition 
on  axial  flow  turbomachine  blade  rows  is  evaluated.  Their  predictive 
equations  and  charts  have  been  presented  as  a  reference  guide  for  deter¬ 
mining  the  transition  region  on  blades.  Additionally,  equations  for  cal¬ 
culating  transitional  boundary  layer  growth  are  presented. 
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2.  GENERAL  TRANSITION  MECHANISM 

Laminar  boundary  layers  formed  on  flat  plates,  curved  surfaces, 
airfoils,  or  compressor  cascades  exhibit  a  strong  tendency  to  break  down 
into  turbulence.  Many  factors,  including  surface  roughness,  freestream 
disturbances,  and  pressure  gradient,  contribute  to  this  breakdown 
process.  However,  the  principal  cause  of  laminar-to-turbulent  transi¬ 
tion  is  the  instability  of  laminar  boundary  layers  to  disturbances  at 
sufficiently  large  Reynolds  number. 

According  to  White  (1) ,  transition  is  principally  a  development  of 
the  boundary  layer  from  instability,  where  infinitesimal  disturbances  grow, 
to  transition,  where  full-scale  turbulence  develops.  However,  White  (1) 
states  emphatically  that  no  theory  of  transition  exists,  and  that  tran¬ 
sition  prediction  relies  on  empirical  and  semi -empirical  correlations 
of  experimental  data.  Morkovin  (2)  reinforces  this  observation  in  his 
uncertainty  principle  for  transition: 

"We  shall  never  know  enough  about  conditions  and  factors 
determining  transition  and  a  basic  uncertainty  will 
remain  for  accurate  quantitative  predictions." 

2.1.  Experimental  Detection  of  Transition 

A  number  of  criteria  are  used  in  boundary  layer  studies  to  locate 
the  transition  region  for  airfoils  experimentally.  According  to 
McDonald  and  Fish  (3),  transition  begins  at  the  point  of  minimum  skin 
friction  coefficient,  C^.  Hall  and  Gibbings  (4)  use  the  minimum  static 
pressure  point  as  measured  with  a  static  pressure  tap.  Schlichting  (5), 
however,  suggests  other  criteria  since  the  minimum  pressure  point  is 


4 


near  transition  only  at  Rec  approximately  10^.  For  larger  values  of 
Reynolds  number,  transition  may  occur  upstream  of  minimum  pressure  (see 
Figure  1).  Also,  at  smaller  values  of  the  lift  coefficient  C  ,  locating 

L 

the  pressure  minimum  may  become  difficult  because  of  "flatter"  pressure 
distribution  along  the  chord.  Schlichting  suggests  that  the  sharp 
decrease  in  the  shape  factor  H,  as  observed  for  flat  plates,  or  a  marked 
increase  in  the  stagnation  pressure  may  give  a  better  indication. 

A  better  indication  is  given  by  the  intermittency  factor  y.  Inter- 
mittency  factor  is  defined  as  the  amount  of  time  the  boundary  layer  is 
turbulent  during  a  given  time  span,  with  y  =  1  for  fully  turbulent  flow. 
It  is  believed  that  the  velocity  distribution  fluctuates  between  a 
laminar  and  turbulent  profile  during  this  time  (5).  Dunham  (6)  charac¬ 
terizes  transition  as  that  region  where  the  intermittency  increases  from 
0.25  to  0.75. 

When  measured  by  hot-wire  anemometers,  the  transition  region  is 
characterized  by  a  series  of  intermittent  turbulent  "bursts."  These 
bursts  take  the  form  of  high  frequency  velocity  fluctuations  in  the 
range  of  1  to  10^  Hz,  with  a  mean  of  about  1000  Hz.  Amplitudes  vary  from 
5  to  20%  of  the  freestream  velocity  (7).  These  turbulent  bursts,  result¬ 
ing  from  the  coalescence  of  turbulent  spots  (8) ,  are  purely  random  in 
nature,  covering  the  entire  frequency  spectrum  without  any  indication  of 
periodicity  (1)  (see  Figure  2). 

Other  indicators  of  transition  presented  by  Schlichting  include 
an  increase  in  the  skin  friction  coefficient,  an  increase  in  the  dynamic 
pressure  measured  at  a  constant  boundary  layer  vertical  location,  and 
development  of  the  fully  developed  turbulent  flow  velocity  profile. 
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Fig.  2 
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.  Hot  wire  measurements  showing  intermittent  turbulent 
fluctuations  from  White  (1). 
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Walker  (8)  used  a  stethoscopic  technique  to  detect  the  acoustic 
emissions  from  the  formation  of  turbulent  spots.  The  formation  of 
turbulent  spots  constituted  the  initial  stage  of  transition  in  his 
boundary  layer  studies  on  compressor  blades.  Walker  also  employed 
a  China-clay  drying  technique. 

2.2.  Qualitative  Description  of  Transition  Initiation 

The  transition  process  generally  results  from  a  buildup  of  dis¬ 
turbances  in  the  boundary  layer.  According  to  Mack  (9),  transition 
results  from  external  disturbances  interacting  with  the  boundary  layer 
either  directly  or  through  instability  waves.  Reshotko  (10)  describes 
transition  as  the  non-linear  response  of  the  laminar  boundary  layer  to 
a  random  forcing  function,  i.e.,  infinitesimal  disturbances.  Roudebush 
and  Lieblein  (11)  state  that  transition  results  from  finite  amplitude 
disturbances  such  as  freestream  turbulence  or  surface  roughness,  or  from 
the  exponential  growth  of  small  instability  waves. 

2.2.1.  White's  Description 

White  (1)  describes  the  mechanism  of  this  process  qualitatively 
and  is  generally  supported  by  many  other  investigators  in  boundary  layer 
theory.  His  six  step  process  consists  of  the  following  steps  as  the 
flow  develops  downstream. 

1.  Initiation  of  unstable  two-dimensional  Tollmien-Schlichting 

(T-S)  waves.  This  is  the  period  of  initial  instability  develop¬ 
ing  from  the  laminar  boundary  layer  as  predicted  from  linear 
stability  theory.  Linear  stability  theory  locates  the  critical 
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Reynolds  number  where  the  T-S  waves  begin  to  grow.  The  T-S 
waves  exhibit  fairly  long  wavelengths,  with  a  minimum  wave¬ 
length  of  66.  These  waves  travel  at  speeds  of  less  than  0.4  U, 
as  predicted  by  linear  stability  theory.  Obremski,  Morkovin, 
and  Landhal  (7)  state  that  this  linear  region  of  the  flow 
ceases  to  exist  when  the  velocity  fluctuation  exceeds  two  per¬ 
cent  of  the  freestream  velocity. 

2.  Development  of  three-dimensional  instability.  The  growth  of 


the  T-S  waves  begins  to  vary  in  the  spanwise  direction.  Walker 
(8)  suggests  that  the  Tollmien-Schlichting  waves  may  possess 


some  slight  three-dimensionality  which  becomes  amplified  by 
secondary  flows  or  turbulence  in  the  freestream.  Next,  accord¬ 
ing  to  Tani  (12),  the  wave  amplitude  develops  a  spanwise  com¬ 
ponent  which  increases  as  the  flow  progresses  downstream.  A 
spanwise  velocity  component,  which  varies  with  time  and  span- 
wise  distance,  also  develops.  This,  in  turn,  initiates  the 
development  of  streamwise  vortices. 

3.  Development  of  a  peak-valley  system.  The  streamwise  vortices 
begin  to  concentrate  locally  and  form  S-shaped  eddies  in  the 
spanwise  direction  (8).  Knapp  and  Roache  (3)  call  this  process 
"thatching."  A  periodic  streamwise  vorticity  system  develops 
into  counterrotating  vortices.  Alternating  velocity  "peaks" 
and  "valleys"  form  in  the  spanwise  direction  with  maximum 
velocity  fluctuations  at  the  peaks.  According  to  Tani,  a  mean 
velocity  defect  is  produced  at  the  peaks,  and  a  velocity  excess 
at  the  valleys. 
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4.  Shear  layer  development  and  vortex  breakdown.  Large  shear 
areas  develop  in  the  boundary  layer  (see  Figure  3)  caused  by  a 
concentration  of  vorticity  at  the  peaks.  An  inflection  point 
develops  instantaneously  and  periodically  in  the  velocity 
profile  of  this  region.  Once  the  shear  layer  forms,  the 
vortices,  which  have  become  stretched  in  an  S-shape  in  the  span- 
wise  direction,  begin  lo  break  down.  According  to  White  (1), 
these  vortices  continue  to  break  down  into  smaller  and  smaller 
vortices  until  all  periodicities  disappear  and  random  fluctua¬ 
tions  set  in.  van  Driest  and  Blumer  (14)  regard  this  step  as 
the  beginning  of  transition. 

5.  Formation  of  turbulent  spots.  When  randomness  sets  in, 
turbulent  spots  begin  to  burst  forth  near  the  peaks  of  the 
vorticity.  These  spots  develop  into  the  familiar  three- 
dimensional  turbulent  spots.  Obremski  and  Fejer  (28)  divide 
this  growth  into  two  modes:  a  "creative"  mode  characterized 
by  accelerating  growth,  and  a  "convective"  mode  characterized 
by  constant  growth.  Such  factors  as  adverse  pressure  gradient 
or  heat  transfer  can  promote  the  creative  growth  mode  and, 
consequently,  speed  up  the  transition  process.  A  number  of 
authors,  including  Walker  (8),  Forest  (15),  and  Dunham  (6), 
regard  this  step  as  the  actual  beginning  of  transition. 

Obremski,  Morkovin  and  Landhal  (7)  have  noted  that  turbulent 
spots  form  when  u’ /U  exceeds  12%.  Along  with  Tani  (12),  they 
note  that  non-linear  regions  act  independently  of  Reynolds  num- 
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6.  Coalescence  of  spots  into  turbulent  flow.  The  flow  between 
the  turbulent  spots  is  still  laminar.  Hot-wire 
in  this  region  will  show  intermittent  turbulence, 
in  the  x  direction  until  complete  turbulent  flow  forms, 
this  region,  the  turbulent  boundary  layer  equations 
applicable.  From  the  turbulent  boundary  layer  viewpoint 
Kline  (43),  energy  transfer  occurs  from  the  mean  flow  to 
boundary  layer  fluctuations.  The  extent  of  turbulence  produc¬ 
tion  should  be  affected  by  the  local  stability  of  the  flow 
(inflection  point  in  the  velocity  profile) .  Turbulence  pro¬ 
duction  is  influenced  by  such  factors  as  pressure  gradient  and 
surface  curvature.  The  entire  process  from  step  1  through 
step  6  is  shown  in  Figure  4. 

Knapp  and  Roache  (13)  add  an  interesting  observation  based  on  flow 
visualization  studies  on  an  axisymmetric  ogive  cylinder.  A  relaminari- 
zation  process  occurs  intermittently  in  the  boundary  layer  after  turbulent 
spot  formation.  The  breakdown  of  the  thatched  vortices  into  smaller 
vortices  causes  an  instantaneous  localized  favorable  pressure  gradient 
of  sufficient  magnitude  to  induce  a  brief  flow  relaminarization.  This 
period  is  of  short  duration,  and  the  vortex  breakdown  process  continues 
immediately  afterwards. 

2.2.2.  Morkovln's  Transition  System 

The  entire  process  described  above  indicates  the  highly  unsteady, 
three-dimensional  character  of  nominally  two-dimensional  flat  plate  flow. 
Flow  conditions  such  as  pressure  gradient,  surface  roughness. 
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4.  Transition  process  on  a  flat  plate  from  White  (1) 


three-dimensionality,  and  heat  transfer  act  upon  various  segments  of 
the  process  to  either  amplify  or  damp  the  disturbances.  Morkovin  (2) 
describes  a  more  generalized  concept  to  include  parametric  effects  of 
the  flow  environment.  His  proposed  system  is  based  on  amplification  of 
the  disturbances  of  linear  stability  theory.  The  diagram  in  Figure  5 
describes  this  system,  illustrating  the  freestream  characteristics  acting 
as  an  oscillator,  providing  input  disturbances  to  an  "operation"  ampli¬ 
fier.  The  flow  environment,  composed  of  the  mean  boundary  layer  proper¬ 
ties  of  pressure  gradient,  surface  roughness,  heat  transfer,  etc., 
determines  the  characteristics  of  the  -~mplif' er.  This  "operation" 
amplifier  produces  three-dimensional  d*  slopment  of  the  disturbances 
and  subsequent  vortex  stretching  in  the  spanwise  direction.  Eventually 
this  leads  to  turbulent  spot  product^  n  and  subsequent  transition. 

Morkovin  also  introduces  the  concept  of  "receptivity"  to  describe 
the  effects  of  the  flow  environment  on  amplification  of  disturbances. 

This  concept  accounts  for  the  scale,  range,  and  frequency  of  disturbances 
similar  to  those  in  the  Tollmien-Schlichting  instability  waves.  The 
disturbance  environment  forces  transition  to  occur.  Unfortunately,  the 
exact  nature  of  this  complex  environment  is  rarely  comprehensible;  conse¬ 
quently,  an  accurate  prediction  of  the  operation  amplifier  is  often 
difficult. 

Reshotko  (10)  describes  receptivity  as  the  response  of  the  laminar 
boundary  layer  to  external  forced  oscillations  of  a  frequency  and  speed 
comparable  to  those  of  the  instability  waves.  The  "signature"  of  the 
external  flow  must  be  known.  Obremski,  Morkovin,  and  Landahl  (7) 


describe  receptivity  as  a  transfer  function  of  the  disturbances  which 
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match  the  T-S  wave  spectra.  They  also  state  that  determination  of  the 
disturbance  environment  is  a  formidable  job. 


The  concept  of  "bypasses"  is  also  proposed,  which  allows  the  input 
and  amplification  of  relatively  high-amplitude  disturbances  such  as 
freestream  turbulence  and  flow  unsteadiness.  The  bypass  mechanism  permits 
turbulent  spot  formation  without  Tollmien-Schlichting  amplification. 

An  example  of  this  is  the  transition  in  pipe  flow,  where  stability 
analysis  does  not  predict  transition.  Another  example  is  two-dimensional 
Poiseuille  flow  in  which  transition  occurs  ahead  of  the  flow  instability 
point. 


2.3.  General  Stability  Theory 

The  appearance  and  growth  of  Tollmien-Schlichting  waves,  as  dis¬ 
cussed  in  step  1  above,  can  be  mathematically  predicted  with  reasonable 
accuracy  by  linear  stability  theory — the  only  viable  theory  of  the 
transition  process.  According  to  Mack  (9,  16),  the  mechanism  by  which 
instability  waves  are  produced  is  unknown.  However,  transition  predic¬ 
tion  and  experimental  verification  depend  on  the  observer's  ability  to 
trace  the  origin  of  the  instability.  Stability  theory  cannot  predict 
the  location  of  transition,  but  it  can  establish  unstable  velocity  pro¬ 
files,  isolate  unstable  frequencies,  and  evaluate  parametric  influences. 

Stability  computations  have  been  conducted  on  Blasius  profiles  by 
Brown  (1959),  on  Falker-Skan  profiles  by  Wazzan,  Okamura,  and  Smith 
(1963),  and  on  plane  Poiseuille  flow  by  Shen  (1954).  The  experiments 
of  Dryden,  Schubauer  and  Skramstad  in  the  early  1940s  at  the  National 
Bureau  of  Standards  verified  predictions  of  stability  theory  for  Blasius  flow. 
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2.3.1.  Linear  Stability  Theory 

The  following  discussion  is  based  on  the  writings  of  Schlichting 
(5),  White  (1),  Mack  (9,  16),  and  Obremski,  Morkovin,  and  Landahl  (7). 
The  cornerstone  of  stability  theory  is  the  Orr-Sonunerf ield  equation, 
which  is  derived  from  the  Navier-Stokes  equation  for  two-dimensional 
flow.  It  assumes  small  disturbance  of  the  mean  boundary  layer  velocity 
U(y),  thus  permitting  a  linear  solution  of  the  Navier-Stokes  equation. 
Next,  the  unsteady  stream  function  is  written  for  the  single  oscillation 


of  the  disturbance. 


i(ax-Bt) 


iKx,y,t)  =  4>(y)e 


where  a  =  2tt/X  represents  the  wave  number,  and  B  is  a  complex  number: 


B  -  Br  +  iB1 


We  now  define  the  complex  quantity 


C  =  -  =  C  +  iC, 
a  r  i 


where  C^.  represents  the  wave  propagation  velocity  and  represents  the 
degree  of  damping  or  amplification  from  the  unsteady  stream  function. 

The  oscillating  components  of  the  boundary  layer  may  be  written  and  sub¬ 
stituted  into  the  Navier-Stokes  equation.  The  result  is  the  fourth 
order  ordinary  differential  equation  for  the  disturbance  amplitude  <J>  (y ) : 


(U  -  C)(<|>"  -  ci2<J>)  -  U"<f>  =  ^  ($""  -  2aV  +  a%  (A) 

6 

where  the  primes  denote  differentiation  with  respect  to  the  dimensionless 
coordinate  y/6.  This  is  an  eigenvalue  problem.  If  a  temporal  stability 


y 
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analysis  is  required  to  determine  the  time-dependent  amplification  of  the 
disturbance,  then  the  four  independent  parameters  of  equation  (4)  are  a, 
Re^,  and  C^,.  If  a  spatial  analysis  is  required  to  calculate  the  dis¬ 
turbance  growth  in  the  x  direction,  then  the  parameters  are  Re^,  a^,  a^, 
and  C.  This  distinction  is  summarized  as  follows. 


Temporal  stability 
a  real 

C  comp lex 

>  0  damped 
<  0  amplification 


Spatial  stability 
a  complex 

C  real 

>  0  damped 
<  0  amplification 


The  eigenvalue  problem  is  solved  by  numerical  integration  of  the  Orr- 
Sommerfeld  equation  which  yields  one  eigenfunction  0(y),  one  complex 
eigenvalue  C  =  +  iC^  for  each  pair  of  a  and  Re^.  Solution  of  this 

equation  yields  phase  velocity,  amplitude  ratio,  overall  growth,  unstable 
frequencies,  and  Re6crit-  The  critical  Reynolds  number  is  the  minimum 
Reynolds  number  on  the  neutral  stability  curve  *  0,  below  which  no  dis¬ 
turbances  of  any  wavelength  will  be  amplified.  This  number  is  used  to 
locate  the  position  of  initiation  of  the  Tollmien-Schlichting  waves. 
Typically,  the  most  unstable  wave  numbers  decrease  as  the  square  root  of 
Re^  for  increasing  Reynolds  number.  The  effect  of  pressure  gradient  can 
be  evaluated  by  determining  the  neutral  stability  curve  and  Re^cr using 
Falkner-Skan  profiles. 

Several  general  theorems  have  been  deduced  from  inviscid  stability 


theory,  which  assumes  large  values  of  Reynolds  number.  This  leads  to  a 
second  order  ordinary  differential  equation  from  the  Orr-Sommerf eld 
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equation.  Two  of  the  theorems  which  have  application  to  viscous  flows 
are 


Theorem  [Tollmien  (1929)] — It  is  sufficient  for  insta¬ 
bility  of  boundary  layer  profiles  that  the  velocity 
profile  have  an  inflection  point. 

Theorem  [Rayleigh  (1880)] — The  phase  velocity  of  an 
amplified  disturbance  must  always  lie  between  the 
minimum  and  maximum  values  of  U(y). 


Adverse  pressure  gradients  in  two-dimensional  flows  typically  induce 
inflection  points  in  viscous  velocity  profiles.  The  speed  of  propaga¬ 
tion  of  amplified  disturbances  (phase  velocities)  is  less  than  0.4  U. 


2.3.2.  Types  of  Flow  Instabilities 

Schlichting  (5)  discusses  flow  instabilities  over  curved  surfaces. 
A  change  in  sign  of  the  expression 


du 


d2u  1 
dy2  R  dy 


(5) 


indicates  instability  on  curved  walls.  R  denotes  local  radius  of  curva¬ 
ture  with  R  >  0  denoting  a  convex  and  R  <  0  denoting  a  concave  wall. 
Centrifugal  forces  on  convex  walls  have  a  stabilizing  effect,  while 
centrifugal  forces  on  concave  walls  have  a  de-stabilizing  effect. 

On  concave  surfaces,  a  primary  instability  in  the  form  of  Gortler 
vortices  develops.  Tani  (12)  indicates  that  these  vortices  cause  a 
spanwise  variation  in  boundary  layer  thickness.  Wortmann  (17),  who  con¬ 
ducted  flow  visualization  tests  with  water,  describes  this  instability 
as  counter-rotating  strearawise  vortices  in  the  spanwise  direction  with 
alternating  inf lection-type  velocity  profiles.  Progressing  downstream. 
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the  Gortler  vortices  oscillate  causing  a  lateral  fluctuation  of  the 
velocity  profile  followed  by  the  formation  of  turbulent  spots.  Dryden 
(18)  suggests  that  this  type  of  instability  dominates  over  Tollmien- 
Schlichting  waves  when  6*/R  >  0.00013. 


2.  A  Laminar  Separation 


Laminar  separation  is  important  in  the  fluid  mechanics  of  airfoils 
and  compressor  blading.  There  is,  in  a  sense,  a  "race"  between  laminar 
separation  and  "natural"  transition  on  the  suction  surfaces.  Roudebush 
and  Lieblein  (11)  point  out  that  transition  makes  axial  flow  compressors 
possible  since  the  resultant  turbulent  boundary  layer  is  less  prone  to 
separation  than  the  laminar  boundary  layer. 

The  initiation  of  laminar  separation  is  independent  of  Reynolds 
number.  It  depends  almost  solely  on  the  pressure  gradient  and  occurs 
when  the  pressure  gradient  parameter  falls  below  the  value  proposed  by 
Thwaites: 


m 


9 


£L  dU 
v  dx 


0.082 


(6) 


This  criterion  is  used  by  Walker  (8),  Seyb  (19),  and  Roberts  (20). 

Dunham  (6),  and  others,  use  -0.09  as  a  criterion.  Citavy  and  Norbury 
(21),  and  Kiock  (22),  report  that  increasing  turbulence  intensity  of  the 
freestream  moves  the  separation  point  downstream.  Their  observations 
are  based  on  cascade  tests  involving  artificially  induced  turbulence. 
Tani  (12)  reports  that  separation  bubble  formation  is  possible  only  when 
Re^  >  500  on  airfoils,  a  claim  verified  by  Owen  and  Klanfer. 
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Laminar  separation  is  undesirable  since  it  causes  high  losses.  It 
is  experimentally  indicated  by  a  sharp  increase  in  momentum  thickness, 

0,  and  a  sharp  decrease  in  shape  factor,  H  (11).  According  to  Roberts 
(20),  separation  is  likely  in  laminar  flow  for  adverse  pressure  gradient, 
smooth  surface,  and  low  Tu  levels.  Five  possible  events  can  occur  in 
relation  to  potential  separation: 

1)  natural  transition  (no  separation) 

2)  short  bubble  formation 

3)  long  bubble  formation 

4)  turbulent  separation 

5)  complete  separation;  no  reattachment 

The  two  significant  events  for  transition  studies  are  the  short 
and  long  bubble  formation.  Roberts  (20)  differentiates  between  the  two 
concepts  when  he  states  that  a  short  bubble  has  a  negligible  effect  on 
the  streamwise  pressure  distribution,  while  a  long  separation  bubble 
affects  the  exterior  freestream  flow  and  the  resultant  pressure  distri¬ 
bution,  including  location  of  the  pressure  peak.  Laminar  separation,  in 
general,  causes  an  initially  thicker  turbulent  boundary  layer  than 
natural  transition. 

A  summary  of  the  laminar  separation  mechanism  will  be  presented, 
based  on  compressor  tests  by  Walker  (8,  23),  cascade  test  data  compiled 
by  Roberts  (20) ,  and  isolated  airfoil  experiments  by  Gault  (24) .  A 
laminar  separation  diagram  based  on  the  work  of  Walker  (8)  is  shown  in 
Figure  6. 

After  laminar  separation  occurs,  a  free  shear  layer  forms  between 
the  top  of  the  boundary  layer  and  the  bubble.  In  the  bubble  itself,  a 
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MEAN  FLOW  DIVIDING  STREAMLINE 

LINE  OF  ZERO  VELOCITY  PARALLEL  TO  WALL 


Fig.  6.  Flow  mechanism  of  a  two-dimensional  laminar  separation 
bubble  from  Walker  (8). 
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small  region  of  reverse  flow  gives  rise  to  a  "dead  air"  region  of  con¬ 
stant  streamwise  pressure.  Dryden  (18)  notes  that  laminar  shear  layers 
are  unstable,  encouraging  transition.  Tani  (12)  suggests  that  this 
occurrence  can  be  approximated  by  the  following  formula: 

l1  =  1006*  (7) 

s 

Horton  (42)  proposes  the  following  criterion  for  airfoils: 

l./eo  *  4  x  104/Reft  (8) 

where  Reg  g  is  the  momentum  thickness  Reynolds  number  at  the  separation 
point.  Tani  and  Roberts  suggest  that  transition  occurs  at  a  particular 
point  while  Walker  (8)  defines  a  finite  length  region.  In  the  transi¬ 
tion  region,  a  turbulent  shear  layer  forms,  a  layer  containing  more 
diffusion  energy  than  the  laminar  shear  layer.  Subsequent  reattachment 
occurs  and  the  amount  of  reverse  flow  in  the  bubble  increases.  Gault 
characterizes  the  turbulent  shear  layer  as  a  region  of  high  vorticity 
which  reattaches  the  flow.  Horton  (42)  determined  that  reattachment 
occurs  under  the  following  condition: 

!£r  =  -°-0059  <9> 

where  R  denotes  the  reattachment  point.  Horton  also  presents  a  semi- 
empirical  equation  for  determining  the  length  of  the  transition  region, 

V 

Pressure  rises  abruptly  at  the  end  of  the  transition  region.  If  a 
long  bubble  forms,  the  turbulent  shear  layer  must  do  more  diffusion  to 
reattach,  causing  lower  pressure  recovery,  reduced  minimum  pressure  on 


< 
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the  suction  surface,  and  reduced  circulation. 

Theoretically,  separation  is  indicated  by  a  zero  wall  shear  stress. 
In  unsteady  flows  in  turbomachinery,  this  condition  may  be  difficult  to 
locate  experimentally. 


2.5.  General  Prediction  Methods 

A  number  of  transition  prediction  criteria  have  been  advanced  based 
on  flat  plate  studies.  They  generally  consist  of  experimental  correla¬ 
tions,  but  some  are  based  on  linear  stability  theory.  A  more  detailed 
treatment  of  these  methods  is  presented  by  Walker  (8),  Morkovin  (2),  and 
Hall  and  Gibbings  (4). 


2.5.1.  The  van  Driest  and  Blumer  Method 

In  the  model  proposed  by  van  Driest  and  Blumer  (14),  a  vorticity 
Reynolds  number  is  defined: 


2  . 

T  =  £_ 
r  V  dy 


(10) 


where  du/dy  is  measured  at  the  point  of  maximum  vorticity,  based  on  the 
observation  that  the  ratio  of  local  inertial  stress  to  local  viscous 
stress  must  reach  a  limiting  value  at  transition.  The  resultant  vorticity 
Reynolds  number  will  reach  a  limiting  value  at  some  value  Re^,  which  must 
be  determined  experimentally.  Effects  of  freestream  turbulence,  surface 
roughness,  and  adverse  pressure  gradient  will  increase  the  value  of  the 
local  vorticity  Reynolds  number  and  reduce  the  transition  Reynolds  number. 
A  curve  fit  was  applied  to  the  data  of  Dryden,  Hall  and  Hislop,  and 


Schubauer  and  Skramstad;  the  correlation  was  then  extended  to  the 
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Pohlhausen  velocity  profiles  to  account  for  the  effect  of  pressure 
gradient : 

(ID 

2.5.2.  The  Crabtree  Method 

Crabtree's  method  (25)  establishes  a  correlation  of  Re„  vs.  m_ 

0tr  9tr 

for  experimental  data  at  low  turbulence,  Tu  <  0.2%.  He  uses  both  flat 
plate  and  isolated  airfoil  data,  and  locates  the  laminar  separation 
region  based  on  Thwaites  criteria.  Crabtree  defines  transition  as  the 
limit  of  existence  of  laminar  flow.  Consequently,  according  to  Hall 
and  Gibbings  (4),  he  is  actually  locating  the  end  of  transition.  To 
use  this  method,  it  is  necessary  to  calculate  the  boundary  layer  growth 
at  each  x  station  and  plot  Reg  vs.  n^.  Transition  is  predicted  at  the 
intersection  of  the  curve  shown  in  Figure  7. 

2.5.3.  The  Granville  Method 

Hall  and  Gibbings  (4)  consider  Granville's  method  very  similar  to 

Crabtree's.  Granville,  however,  plots  (RQ  -  Rrt  . _)  vs.  m„,  where  mQ 

r  0tr  0crit  0  0 

is  the  average  pressure  gradient  parameter  within  the  instability  region. 
In  other  words,  Granville's  method  is  basically  a  replotting  of  Crabtree's 
data.  Granville  studied  the  effect  of  freestream  turbulence  and  showed 
that  Reg  ^  approaches  a  value  of  180  as  Tu  exceeds  5%. 

2.5.4.  The  e11  Method 

This  method  is  based  on  a  strict  extension  of  stability  theory 
computations.  The  ratio  of  the  disturbance  amplitude  to  the  initial 


=  1  -  0.0485  X  +  3.36  Re 


6,tlu 


til 


7.  Transition  prediction  data  correlation  from  Crabtree  (25). 
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value  is  computed.  Van  Ingen  (29)  computed  these  ratios  from  stability 
diagrams  and  compared  them  to  experimental  results  on  isolated  airfoils. 
The  amplification  exponent  is,  by  definition, 

n  »  In  ~  (12) 

o 

where  Aq  is  the  amplitude  of  the  initial  disturbance.  Van  Ingen  main¬ 
tained  that  stability  theory  is  good  for  predicting  the  early  stages  of 
transition.  Basing  his  method  on  airfoil  experiments,  he  proposed  a 
value  of  n  *  7  or  8.  Smith  and  Gamberoni  recommended  n  =  9.  Van  Ingen 
did  note  that  the  method  becomes  inaccurate  in  adverse  pressure  gradients 
near  the  point  of  laminar  separation. 

All  of  the  methods  mentioned  thus  far  require  computing  the  boundary 
layer  variables  6,  0,  6*,  Re^  at  each  station  and  comparing  them  to  the 
critical  values  or  transitional  values  at  that  station. 

2.5.5.  The  Hall  and  Gibbings  Method 

Hall  and  Gibbings  (4)  present  a  general  criterion  for  locating 
transition  based  on  the  data  correlation  of  flat  plate  and  concave 
surface  studies.  Their  correlation  (see  Figure  8)  is  useful  since  the 
data  represent  different  values  of  pressure  gradient.  Compressor  blade, 
cascade,  and  isolated  airfoil  data  involve  variable  pressure  gradient 
in  the  streamwise  direction. 

Hall  and  Gibbings  define  the  transition  location  as  the  point  of 
minimum  pressure.  As  previously  mentioned,  this  is  not  a  very  accurate 
method  because  of  the  Reynolds  number  effect. 


1  -1— *  t  a  -  -*-■>  ■’  « 
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3.  PARAMETRIC  EFFECTS  ON  TRANSITION 

In  this  section,  parameters  which  affect  the  transition  process 
will  be  examined  with  emphasis  placed  on  those  variables  having  particu¬ 
lar  application  to  flow  thru  axial  flow  compressors,  cascades,  and 
isolated  airfoils.  A  review  of  the  pertinent  literature  revealed  that 
while  extensive  tests  have  been  conducted  on  individual  effects,  combined 
effects  have  not  been  adequately  investigated. 

The  major  influences  on  transition  are  the  pressure  gradient, 
turbulence  intensity,  and  surface  roughness,  though  surface  curvature, 
heat  transfer,  and  acoustic  disturbances  also  have  some  effect.  The 
relative  importance  of  each  parameter  depends  on  the  particular  flow 
situation. 


3.1.  The  Effect  of  Pressure  Gradient 

The  freestream  pressure  gradient  has  the  greatest  effect  on  transi¬ 
tion.  Walker  (8)  points  out  that  pressure  gradient  effects  account  for 
about  one-half  of  all  combined  effects  for  two  reasons.  First,  the 
instability  point  Re^^  is  influenced  strongly  by  pressure  gradient; 
second,  the  distance  between  transition  and  instability  (Re5*tr  “  Re5*i^ 
is  affected.  In  this  region,  the  amplification  factors  are  altered  by 
pressure  gradient. 

The  external  pressure  gradient  for  airfoils  and  airfoil  cascades 
can  be  calculated  by  several  techniques  using  iterative  inviscid/boundary 
layer  computation.  These  techniques  are  summarized  in  a  review  by 
Hansen  (27). 


PRBCKD1N0  PiOft  BUMUNOT  FILMS 


30 


The  effect  of  pressure  gradient  is  stronger  on  the  instability 
Reynolds  number  than  on  the  transition  Reynolds  number.  White  (1) 
obtained  the  following  approximations  for  Pohlhausen  profiles: 

ReS*crit  ’  520  »°'6X  (13> 

Re6*cr  ^  2,00  e°'°8)  <14> 

Reshotko  (10)  points  out  that  a  favorable  gradient  has  a  stronger 
effect  on  instability  than  an  adverse  pressure  gradient.  This  phenomenon 
is  illustrated  for  Falkner-Skan  laminar  profiles  in  Figure  9.  The  over¬ 
all  effect  of  pressure  gradient  on  transition  is  shown  in  Figure  10. 

adverse  pressure  gradients  tend  to  act  on  the  step  3  mechanism  of 
transition.  The  S-shaped  spanwise  vortices  show  less  tendency  to 
"thatch"  and  are  more  rounded,  indicating  less  vortex  stretching  and 
more  compression.  The  compression  of  these  vortices  causes  more  rapid 
breakdown  (13)  and  promotes  the  "creative"  growth  mode  of  turbulent  spot 
production  (28). 


3.2.  The  Effect  of  Freestream  Disturbances 

The  effects  of  freestream  turbulence  are  perhaps  the  most  contro¬ 
versial  aspect  of  boundary  layer  transition.  Conflicting  ideas  have 
arisen  regarding  their  effect  on  stability  and  on  final  transition. 
Freestream  turbulence  in  axial  flow  turbomachinery  consists  of  two  com¬ 
ponents:  random  fluctuations  of  the  mean  flow  and  flow  unsteadiness  due 
to  the  stator  and  rotor  wakes.  It  is  worthwhile  to  consider  the  effect 


•  c 


of  unsteadiness  effects  separately. 


3.  PARAMETRIC  EFFECTS  ON  TRANSITION 


In  this  section,  parameters  which  affect  the  transition  process 
will  be  examined  with  emphasis  placed  on  those  variables  having  particu¬ 
lar  application  to  flow  thru  axial  flow  compressors,  cascades,  and 
isolated  airfoils.  A  review  of  the  pertinent  literature  revealed  that 
while  extensive  tests  have  been  conducted  on  individual  effects,  combined 
effects  have  not  been  adequately  investigated. 

The  major  influences  on  transition  are  the  pressure  gradient, 
turbulence  intensity,  and  surface  roughness,  though  surface  curvature, 
heat  transfer,  and  acoustic  disturbances  also  have  some  effect.  The 
relative  importance  of  each  parameter  depends  on  the  particular  flow 
situation. 


3.1.  The  Effect  of  Pressure  Gradient 

The  freestream  pressure  gradient  has  the  greatest  effect  on  transi¬ 
tion.  Walker  (8)  points  out  that  pressure  gradient  effects  account  for 
about  one-half  of  all  combined  effects  for  two  reasons.  First,  the 
instability  point  Re^^  is  influenced  strongly  by  pressure  gradient; 
second,  the  distance  between  transition  and  instability  ( Re <5 * t r  “ 
is  affected.  In  this  region,  the  amplification  factors  are  altered  by 
pressure  gradient. 

The  external  pressure  gradient  for  airfoils  and  airfoil  cascades 


can  be  calculated  by  several  techniques  using  iterative  inviscid/boundary 
layer  computation.  These  techniques  are  summarized  in  a  review  by 
Hansen  (27). 
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Neutral  stability  diagrams  for  Falkner-Skan  flows  showing  the 
effect  of  pressure  gradient  from  Reshotko  (10).  Parallel  flow 
indicates  a  streamwise  component  of  the  mean  flow  velocity. 
Non-parallel  indicates  existence  of  a  vertical  component. 
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Fig.  10.  Transition  data  correlation  showing  effect  of  pressure 
gradient  from  Hall  and  Gibbings  (4).  \q  denotes  the 
average  pressure  gradient  parameter  over  the  instability 
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3.2.1.  The  Effect  of  Flow  Unsteadiness 

Data  on  periodic  unsteadiness  comes  principally  from  Obremski  and 
Fejer  (28) ,  who  conducted  flat  plate  wind  tunnel  tests  with  oscillating 
shutter  blades.  They  defined  transition  as  the  location  where  the  inter- 
mittency  factor  y  >  0.  Frequency  of  the  oscillations  ranged  from  5  to 
60  Hz,  with  Tu  <  .2%.  They  defined  a  "non-steady"  Reynolds  number 


Re 


Ns 


L*AU 

2itv 


(15) 


where  L*  is  a  characteristic  length  parameter.  If  this  value  exceeded 
27,000,  the  transition  process  was  periodic;  if  less  than  25,000,  the 
transition  was  "aperiodic."  Generally,  unsteadiness  had  little  contri¬ 
bution  in  favorable  pressure  gradients  if  AU/U  <  0.20.  For  adverse 
pressure  gradients,  ReNs  >  6,000  indicated  periodicity. 

Downstream,  turbulent  bursts  appear  in  the  troughs  of  the  unsteady, 
instantaneous  velocity  trace.  The  flow  unsteadiness  appears  to  promote 
the  "creative"  mode  of  turbulent  spot  growth  with  subsequent  reduction 
in  transition  lengths.  Steady  flow  favors  the  "convective"  mode  of  spot 
development. 

Dzung  (29)  also  points  out  that  unsteady  flow  generates  a  "starting 
vortex"  on  each  cycle  because  of  the  time-varying  circulation.  This 
vortex  is  convected  downstream  to  the  next  set  of  blade  rows.  Walker  (8) 
reinforces  these  observations  in  his  compressor  studies  and  adds  that 
when  Re^g  exceeds  the  critical  value,  the  spots  appear  in  ribbon-like 
form  in  the  spanwise  direction. 

Obremski  and  Fejer  also  noticed  a  "calm"  period  after  the  passing 


of  turbulent  spots  similar  to  the  calm  period  c  served  by  Roache  and  Knapp  (13) . 
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3.2.2.  The  Effect  of  Overall  Turbulence  Intensity 

The  turbulence  intensity  parameter  Tu  consists  of  two  parts:  a 
periodic  flow  unsteadiness  term  and  a  random  velocity  fluctuation  term. 
According  to  the  single-stage  compressor  test  results  of  Evans  (30) , 
these  quantities  are  of  equal  magnitude  at  design  speed.  The  periodic 
unsteadiness  frequencies  are  equal  to  the  blade  passing  frequencies. 

Random  fluctuations  have  frequencies  in  the  range  of  20  Hz  to  20  KHz  (22). 

Using  compressor  data,  Kiock  (22)  reports  that  turbulence  increases 
rapidly  from  1%  to  10%  through  a  three-stage  compressor.  Schlichting 
and  Das  (31)  quote  other  sources  showing  an  increase  of  6  to  8%  through 
multi-stage  compressors.  Laminar  separation  bubbles  tend  to  increase 
the  downstream  freestream  turbulence.  Citavy  and  Norbury  have  observed 
a  small  decrease  in  the  pressure  coefficient  distribution  across  the 
cascade  chord  as  turbulence  increases.  This  change  may  be  due  to  changes 
in  the  laminar  separation  bubble  or  to  the  axial  velocity  ratio  which 
was  also  varied. 

Using  cascade  test  data,  Schlichting  and  Das  (31)  observed  that  as 
turbulence  increased  from  0.5%  to  2%  there  was  an  increase  in  boundary 
layer  momentum  thickness  across  the  cascade  chord.  They  also  observed 
a  decrease  in  the  shape  factor.  Little  change  was  observed  above  2% 
turbulence,  indicating  formation  of  a  turbulent  boundary  layer.  This 
effect  was  also  observed  by  Schaffler  (32). 

Schlichting  (5)  presents  a  strong  case  for  turbulence  effects  on 
transition  in  a  zero  pressure  gradient.  Figure  11  shows  that  turbulence 
has  little  effect  on  duration.  Though  stability  aspects  of  turbulence 
are  not  fully  understood,  Obremski  (7)  feels  that  freestream  turbulence 
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Effect  of  freestream  turbulence  on  the  transition- 
instability  location  in  a  zero  pressure  gradient,  from 
Schlichting  (5). 
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may  contain  frequencies  that  are  in  the  range  of  the  Tollmien-Schlichting 
selective  amplifier,  while  Morkovin  (2),  in  his  receptivity  theory, 
allows  turbulent  disturbances  of  sufficient  magnitude  to  "bypass"  the 
Tollmien-Schlichting  mechanisms  altogether. 

Rogler  (33)  presents  the  most  cogent  description  of  turbulence 
effects  on  stability.  Freestream  disturbances  exist  as  vortex  sheets, 
as  rings  and  patches,  or  as  wakes  which  convect  into  the  boundary  layer, 
an  area  of  concentrated  vorticity.  A  spatial  analysis  is  conducted  by 
writing  a  two-dimensional  vorticity  equation  comprising  convection,  pro¬ 
duction,  and  diffusion  terms  for  the  boundary  layer.  The  ingested  free¬ 
stream  vorticity  then  induces  three-dimensional  velocities  which  generate 
additional  boundary  layer  vorticity.  A  readjustment  of  the  vorticity 
occurs,  varying  with  time  or  distance.  This  disruptive  activity  induces 
forced  oscillations  with  wavelengths  in  the  range  of  the  Tollmien- 
Schlichting  amplifier  followed  by  turbulent  spot  formation. 

On  a  macroscopic  level,  Walker  (8,  23)  presents  an  entirely  dif¬ 
ferent  view;  he  states  that  turbulence  influences  transition  indirectly 
by  directly  modifying  the  pressure  distribution.  This  assertion  is 
supported  by  the  cascade  tests  of  Citavy  and  Norbury  (21)  and  Kiock  (22) . 
Walker  states  that  large  scale  freestream  disturbances  cause  instantaneous 
changes  in  flow  incidence  which  move  the  instantaneous  suction  peak  and 
subsequent  transition  location.  Walker  then  states  that  freestream 
turbulence  has  an  effect  on  transition  only  if  turbulence  scale  alters 
the  pressure  distribution  and  incidence,  changes  the  initial  amplifica¬ 
tion  of  neutral  disturbances,  or  changes  the  amplification  rate  of  all 
disturbances  wavelengths.  If  none  of  these  three  conditions  are  met. 
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turbulence  intensity  is  not  expected  to  move  up  the  transition  location 
by  more  than  20%  of  xt  -  x^.  Proponents  of  stability  theory  admit  that 
the  effect  of  turbulence  intensity  on  instability  is  largely  unknown 
(7);  consequently.  Walker's  observations  have  some  basis. 

Horlock  and  Evans  (34)  attempted  a  theoretical  treatment  of  turbu¬ 
lence  intensity  through  use  of  the  von  Karman  momentum  integral  equation. 
They  suggest  that  additional  terms  be  included  in  the  equation  and  that 
the  shear  stress  does  not  approach  zero  at  the  edge  of  the  boundary 
layer.  They  propose  the  following  relation: 

T. 

=  4  Tu  (16) 

wall 


3.3.  The  Effect  of  Surface  Roughness 

Much  of  the  work  on  effects  of  surface  roughness  (12)  has  been 
conducted  on  two-dimensional  tripwires  or  three-dimensional  isolated 
elements,  but  it  is  not  entirely  clear  which  type  of  element  applies  to 
two-dimensional  boundary  layers. 

First  of  all,  Schaffler  (32)  states  that  it  is  difficult  to  charac¬ 
terize  roughness  because  of  the  individual  shape  of  peaks,  the  roughness 
peaks  per  unit  of  length,  and  the  uniformity  of  roughness  heights.  The 
arithmetic  average  (AA)  alone  does  not  give  a  good  description  of  the 
hydrodynamic  characteristics.  Schaffler  defines  the  height  parameter  k 
as  the  difference  between  the  10  highest  peaks  and  10  deepest  valleys  in 
100  pm  of  length.  Thus,  his  roughness  criterion  is 
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Re,  100 

k,crit  v 


for  hydrodynamically  smooth  surfaces  where  W  is  the  free  stream  velocity 
relative  to  the  cascade. 


For  a  single  two-dimensional  roughness  element,  Schlichting  recom¬ 


mends 


Ukk 

Re,  =  -i-  =  20 

k,crit  v 


For  an  isolated  three-dimensional  element  on  two-dimensional  flow, 
Tani  (12)  recommends 


\k 

Rek,crit  "  -T-  =  600 


Tani  also  notes  that  pressure  gradient  has  no  effect  on  the  value  of 

Re,  . _  in  three-dimensional  flow, 
k.crit 

For  distributed  roughness,  which  is  most  applicable  to  actual  turbo¬ 
machinery  flows,  Schlichting  (5)  recommends 


Re,  . - r  =  120 

kcrit  V 


Below  this  number,  roughness  has  no  effect  on  transition.  In  Figure  12 
Schlichting* s  graph  shows  the  effect  of  surface  roughness  and  pressure 
gradient. 

McDonald  and  Fish  (3)  present  a  convenient  graph  for  distributed 
roughness  based  on  the  parameter  k/6*  (see  Figure  13). 

Referring  to  stability  theory,  Morkovin  (2)  states  that  roughness 
does  not  generate  unsteady  disturbances,  but  rather  acts  as  a  modifier 
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MEASUREMENTS  OF  FEINDT  (REF.  32),  SANDPAPER  ROUGHNESS,  FOR  VARIOUS 
STREAMWISE  PRESSURE  GRADIENTS 
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Fig.  13.  Effect  of  distributed  surface  roughness  on  transition 
Reynolds  number  from  McDonald  and  Fish  (3). 
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to  the  instability  waves.  Dryden  (18)  takes  the  opposite  view  and  argues 
that  roughness  elements  generate  disturbances  through  eddy  production. 


3.4.  The  Effect  of  Heat  Transfer  and  Compressibility 

In  incompressible  flow,  heat  transfer  to  the  fluid  causes  ^ecr^t  to 
decrease  because  of  the  dependence  of  viscosity  on  temperature;  conversely, 
heat  transfer  from  the  fluid  has  a  stabilizing  effect.  Schlichting  (5) 
demonstrates  mathematically,  using  the  Navier-Stokes  equation,  how  heat 
transfer  to  the  fluid  generates  an  inflection  point  in  the  velocity  pro¬ 
file. 

In  compressible  flow,  a  thermal  boundary  layer  develops.  Since  the 
velocity  at  the  wall  is  zero,  the  wall  temperature  approaches  the  stag¬ 
nation  temperature,  as  indicated  by  the  recovery  factor.  For  subsonic 
flow,  this  resultant  heat  transfer  effect  is  slight  up  to  a  Mach  number 
of  0.9. 

At  higher  Mach  numbers,  however,  local  shock  waves  may  be  produced, 
generating  a  strong  adverse  pressure  gradient  and  inducing  transition 
(32).  The  study  of  transition  in  supersonic  flows  is  of  great  interest 
in  aeronautical  engineering  and  a  wealth  of  experimental  data  is  avail¬ 
able,  though  these  applications  are  beyond  the  scope  of  this  paper. 

3.5.  The  Effect  of  Streamwise  Surface  Curvature 

Centripetal  forces  generated  on  concave  walls,  such  as  pressure 
surfaces  on  airfoils,  give  rise  to  a  primary  flow  instability  called 
Gortler  vortices  which  are  characterized  by  alternating  streamwise 
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vortices  of  opposite  rotation  and  alternating  inflection-type  two-dimen¬ 
sional  velocity  profiles  (secondary  instability)  (see  Figure  14) .  In 
extensive  flow  visualization  studies  of  these  effects  done  by  Wortmann 
(17),  the  vortices  show  an  oscillating  motion  and  the  velocity  peaks  and 
valleys  show  a  lateral  movement.  This  instability  destroys  the  symmetry 
of  the  vortices,  and  a  fluctuation  of  the  velocity  profile  occurs  in  all 
three  planes,  followed  by  turbulent  spot  formation. 

Dryden  (18)  indicates  that  Gortler  vortices  become  dominant  if  6*/r 
>  0.00013.  Tani  (12)  locates  the  point  of  instability  through  a  Gortler 
number,  and  states  the  following  stability  limit: 

G  =  Rg*(fi*/r)1^2  >  1.2  (21) 

Tani  also  presents  a  diagram  (see  Figure  15)  to  determine  the  transition 
Reynolds  number. 

On  convex  surfaces,  centrifugal  forces  arising  from  the  suction 
surface  curvature  exert  a  stabilizing  influence  on  transition  (5).  If, 
however,  6/r  is  much  less  than  1,  this  influence  is  very  small. 
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Fig.  15.  Effect  of  surface  curvature  on  transition  in  a  zero  pressure 
gradient  from  Tani  (12) . 
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Fig.  16.  Variation  of  the  transition  length  in  a  laminar  separation 
bubble  with  turbulence  factor  from  Roberts  (20). 
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4.  TRANSITION  IN  TURBOMACHINE  BLADE  ROWS 

The  discussion  will  now  focus  on  boundary  layer  flows  through  axial 
flow  turbomachine  blade  rows  with  emphasis  on  compressor  blades,  particu¬ 
larly  the  suction  surface.  Boundary  layers  form  over  the  hub  and  outer 
walls  as  well,  but  these  phenomena  will  not  be  discussed  here.  The 
importance  of  accurately  determining  the  location  of  transition  will  be 
considered,  followed  by  a  discussion  of  the  boundary  layer  phenomena 
which  occur  over  the  blades  and  several  prediction  techniques. 

4.1.  The  Importance  of  Accurately  Locating  Transition 

Location  of  transition  is  important  in  the  study  of  losses  in 
turbomachinery.  According  to  Koch  and  Smith  (35),  loss  sources  include 
blade  profile  losses,  end  wall  boundary  layers,  shock  losses,  and  drag 
losses  from  part  span  shrouds.  Blade  profile  losses  are  a  boundary 
layer  phenomenon  which  generates  a  stagnation  pressure,  P  ,  loss  through 
the  blades  and  creates  wakes  (vorticity) .  Both  of  these  effects  generate 
entropy  with  a  resultant  loss  of  thermal  efficiency. 

Lieblein  (36)  observes  that  losses  arise  from  the  growth  of  the 
boundary  layer.  Also,  Pq  losses  can  be  correlated  to  the  boundary  layer 
momentum  thickness  and  shape  factor  at  the  trailing  edge.  Since  the 
boundary  layer  growth  characteristics  are  quite  different  for  laminar, 
transitional,  and  turbulent  boundary  layers,  analytic  loss  predictions 
depend  on  an  accurate  calculation  of  the  location  and  extent  of  each 
region.  Roudebush  and  Lieblein  (11)  estimate  that  the  calculation  of 
the  drag  coefficient  may  be  in  error  by  10%  because  of  an  error  in  location 
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of  transition  by  20%  of  chord  at  Re  =  10^. 

4.2.  Parametric  Effects  in  Turbomachinery  Flows 

Typical  boundary  layer  effects  observed  in  axial  flow  compressors 
and  the  resultant  effects  on  transition  are  summarized  below.  This 
summary  includes  the  effects  of  freestream  turbulence,  surface  roughness, 
and  laminar  separation. 

4.2.1.  Reynolds  Number  Effects 

Schlichting  and  Das  (31)  summarize  the  effect  of  chord  Reynolds 

number  Rec  as  follows:  for  Rec  <  5  x  10^,  transition  is  significant  on 

the  suction  surface;  for  Rec  <  2  x  10^,  laminar  separation  begins;  for 

Rec  <  1  x  lo\  complete  laminar  separation  is  expected  to  occur.  These 

results  are  for  low  turbulence  levels  Tu  <  2%  based  on  cascade  tests. 

In  compressor  tests  where  Rec  =  1.5  x  10^,  Walker  (8)  observed  that 

the  transition  region  occupied  between  15%  to  45%  of  total  chord  length. 

The  start  of  transition  was  independent  of  Re^ ,  but  the  instability 

region  ReQ.  -  ReQ .  decreased  as  Re  decreased.  Turbulence  intensities 
Ut  t)l  c 

were  in  the  range  of  2%  to  6%. 

Pollard  and  Gostelow  (37)  observed  laminar  separation  on  cascade 

tests  for  Re  =  1.1  x  10^  with  low  turbulence  levels  and  noted  that 
c 

decreasing  Rec  caused  an  increase  in  bubble  length. 

Schaffler  (32)  observed  that  turbulence  intensity  had  no  effect  on 

4 

the  length  of  the  separation  bubble  in  cascade  tests  for  Rec  =  5  x  10  . 

4.2.2.  The  Effect  of  Freestream  Turbulence  and  Unsteadiness 
Turbulence  increases  in  intensity  as  the  flow  progresses  through 


49 


each  stage  of  a  compressor.  Dunham  (6)  presents  an  equation  to  approxi¬ 
mate  the  increase  in  turbulence  based  on  the  contraction  ratio.  Turbu¬ 
lence  levels  at  high  flight  altitudes  are  generally  nil  and  may  cause 
laminar  separation  in  the  first  stage  of  compressors. 

Evans  (30,  38)  measured  the  turbulence  and  unsteadiness  in  single- 
stage  compressors.  His  measurements  indicate  a  turbulence  intensity 
(including  random  fluctuation  and  unsteadiness)  of  2.5%  to  7%,  with 
decreasing  turbulence  for  increasing  flow  coefficient.  He  notes  that 
the  velocity  drops  sharply  on  the  pressure  side,  but  rises  more  slowly 
on  the  suction  surface,  indicating  a  thicker  boundary  layer.  He  also 
notes  that  the  tangential  fluctuating  velocity  reaches  a  maximum  in 
phase  with  the  passing  wake.  Velocity  defect  in  the  wake  is  approxi¬ 
mately  20%.  This  is  similar  to  the  findings  of  Roudebush  and  Lieblein. 
Evans  also  estimates  that  the  measured  velocity  profiles  fluctuate 
periodically  at  transition  between  a  fully  laminar  and  fully  turbulent 
profile.  The  largest  unsteady  profile  fluctuation  occurred  at  y/6  =  0.3. 

Walker  indicates  that  the  main  effect  of  periodic  unsteadiness 
caused  by  the  passing  of  wakes  is  "smearing  out"  of  the  transition 
region  (30). 

4.2.3.  Surface  Roughness 

In  three-stage  compressor  tests  with  NACA  65  blades,  Bammert  and 
Woelk  (39)  observed  that  surface  roughness  distributed  on  the  blade 
surface  beyond  a  60  grain  size  caused  a  decrease  in  efficiency  of  2  to 
6%.  Schaffler  reports  a  decrease  of  3%.  Neither  author,  however,  con¬ 
ducted  hot-wire  measurements;  consequently,  no  velocity  profiles  were 
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determined • 

4.2.4.  Laminar  Separation 

Laminar  separation  causes  an  efficiency  loss  in  compressors.  Unfor¬ 
tunately,  a  theoretical  treatment  of  separation  is  not  presently  possible 
The  following  results  are  based  on  the  NACA  65  cascade  tests  of  Roberts 
(20)  and  Kiock  (22),  and  the  compressor  tests  of  Walker  (8,  23).  All 
tests  employ  Thwaite's  separation  criteria. 

Roberts  defines  a  "bursting"  Reynolds  number  Rec(g)  where  the  loss 
coefficient  increases  rapidly  because  of  a  long  separation  bubble. 
Typically,  the  separation  point  is  independent  of  chord  Reynolds  number 
but  dependent  on  angle  of  attack.  Roberts'  correlation  is  based  on  a 
turbulence  factor  CTFI  which  accounts  for  the  scale  of  the  turbulence. 

The  bursting  Reynolds  number  Re^^  =  U^c/v  is  determined  from  the  fol¬ 
lowing  correlations: 

Rec(B)  =  7~:YTfF)^  x  1C)5  +  10’000  (22> 

where  D  is  the  diffusion  factor  related  to  upstream  and  downstream 
velocities  as  defined  by  Roberts  (2).  Roberts  also  establishes  a  cor¬ 
relation  to  determine  the  location  of  transition  in  the  bubble  from  the 
separation  point  (see  Figure  16). 

Kiock  conducted  cascade  tests  using  an  oscillating  grid  to  simulate 
unsteadiness,  varying  the  turbulence  intensity  from  0.3%  to  8%.  He 
concluded  that  at  constant  Rec  turbulence  intensity  reduced  bubble  length 
Increasing  the  chord  Reynolds  number  at  constant  Tu  reduces  bubble 
length.  Laminar  separation  was  observed  for  chord  Reynolds  number  of 


r 
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90,000  <  Re  <  1.6  x  105. 
c 

Walker  uses  Horton's  (42)  criteria  to  determine  the  following  tran¬ 
sition  location: 
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(23) 


where  the  mean  shape  factor  across  the  bubble  is  less  than  4.20.  He 
proposes  a  different  equation  based  strictly  on  a  curve  fit  if  the  shape 
factor  is  greater  than  4.20.  The  turbulent  reattachment  point  criteria 
is  based  on 


0_  dU 
U  dx 


-0.0082 


(24) 


Generally,  laminar  separation  occurred  on  the  British  C4  section  profiles 
below  Rec  <  2  x  10“* .  The  separation  point  moves  downstream  and  increases 
in  bubble  length  as  the  angle  of  incidence  decreases,  supporting  Roberts’ 
observation.  Walker  also  recommends  using  a  linear,  rather  than  a 
parabolic,  curve  fit  to  the  U(x)  data  to  locate  the  point  of  separation. 

4.2.5.  Turbine  Blading 

Turbine  blade  rows  present  a  different  transition  problem  than  com¬ 
pressor  blades  because  of  pressure  gradient  and  heat  transfer  effects. 
Since  turbine  blades  are  often  cooled  internally,  heat  transfer  from  the 
fluid  tends  to  stabilize  the  flow.  Seyb  (19)  notes  that  pressure  dis¬ 
tributions  are  such  that  separation  does  not  occur  unless  induced  by 
shocks.  Bay  ley  (40)  notes  that  for  2  x  10“*  <  Rec  <  2  x  10^  the  pressure 
surface  stays  laminar,  while  suction  surface  transition  occurs  in  the 
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range  of  50%  to  70%  of  chord.  Bayley  asserts  that  turbulence  levels 
for  turbines  are  higher,  usually  around  20%. 

4.3.  Prediction  Methods 

All  of  the  following  methods  for  predicting  transition  are  based  on 
either  direct  experimental  observations  or  compressor/cascade  data  cor¬ 
relations  methods  using  Thwaites'  method  for  calculation  of  the  laminar 
boundary  layer  thickness,  Martensen’s  method  for  calculating  laminar 
velocity  profiles,  and  Truckenbrodt' s  method  for  turbulent  boundary 
layers. 

4.3.1.  General  Methods 

Experimental  techniques  exist  which  give  estimates  of  the  transition 
region,  usually  based  on  hot-wire  anemometer  measurements  in  cascades 
and  compressors.  Transition  causes  a  sharp  decrease  in  the  shape  factor 
H,  which  indicates  formation  of  a  turbulent  velocity  profile.  Alterna¬ 
tively,  transition  can  be  located  by  a  sharp  increase  in  the  skin  friction 
coefficient.  If  a  pitot  tube  is  placed  in  the  boundary  layer,  an 
increase  in  dynamic  pressure  can  indicate  the  beginning  of  the  transition 
region.  Some  investigators  use  the  minimum  suction  pressure  point. 
However,  as  previously  discussed  and  as  illustrated  in  Figure  1,  this 
is  only  valid  at  high  Rr.c  numbers.  A  hot-wire  anemometer  may  be  used  to 
indicate  the  appearance  of  turbulent  bursts,  and  hence  determine  the 
intermittency  factor.  However,  these  bursts  begin  to  appear  in  the 
shear  layer  of  the  boundary  layer,  so  the  investigator  should  estimate 
this  layer  using  the  suggestions  of  Tani  (12)  and  Walker  (8). 
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All  of  the  experimental  techniques  mentioned  above  require  that 
the  observer  understand  the  transition  mechanism  in  order  to  adequately 
analyze  the  data. 

4.3.2.  Dunham's  Correlation 

Dunham's  (6)  estimations  are  based  on  a  curve  fit  of  turbine  blade 

cascade  and  compressor  test  data  from  other  sources.  Parameters  are 

2 

pressure  gradient  m0  =  0  /v(dU/dx)  and  an  average  turbulence  intensity 
between  leading  edge  and  trailing  edge  of  the  blade.  Dunham  safely 
defines  transition  as  the  region  where  the  intermittency  varies  from 
0.25  to  0.75. 

His  correlation  equation  is 

Re0  tr  =  (0.27+0.73  exp (- . 8Tu) )  (550+680  (l+-Tu-21m0))  (25) 

where  Tu  represents  the  average  turbulence  intensity.  Good  correlation 
was  established  with  Evan's  compressor  data,  but  less  success  was  obtained 
with  turbine  data. 

The  author  of  this  paper  recommends  using  this  equation  with  caution 
since  only  one  source  for  compressor  data  was  used  and  a  compressor  with 
a  different  blade  profile  may  yield  different  results.  Dunham's  turbine 
data  are  more  comprehensive  since  several  sources  were  used. 

Dunham  also  suggests  a  correlation  for  predicting  the  transition 
length  in  a  laminar  separation  bubble,  based  on  Walker's  data: 

R  -  =  5000  +  exp (10. 463  -  1.54  Tu0,629)  (26) 

il 

Unfortunately,  Dunham  again  used  only  one  source.  This  author  assumes 
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that  the  compressors  used  by  Walker  and  Evans  were  equipped  with  British 
C4  section  blades;  it  should  be  noted  that  other  profiles  may  have  dif¬ 
ferent  characteristics.  Dunham's  definition  of  transition,  though, 
agrees  well  with  White’s  (1). 


4.3.3.  Seyb's  Correlation 

Seyb  (19)  takes  a  similar  approach  and  arrives  at  a  correlation  in 
the  form  of  a  chart  (see  Figure  17).  This  chart  is  very  similar  to 
Dunham's,  but  gives  no  reference  to  the  sources  of  the  data.  Also, 
Seyb's  chart  does  not  extend  to  the  higher  turbulence  levels  found  in 
turbomachinery.  The  equation  for  this  correlation  is 
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(27) 


for  turbulence  intensity  less  than  4%. 

4,3.4.  Forest's  Method 

Forest  uses  a  turbulence  model  approach  and  accounts  for  accelera¬ 
tion  forces,  curvature,  turbulence,  and  pressure  gradient.  Forest's 
principal  contribution  is  the  development  of  a  scheme  to  calculate  the 
boundary  layer  development  in  the  transition  region,  the  velocity  pro¬ 
files  in  that  region,  and  the  end  of  transition.  He  uses  Seyb's  criterion 
to  predict  the  beginning  of  transition. 

An  effective  viscosity  is  used  to  account  for  increasing  turbulence 
in  the  transition  region,  utilizing  an  intermittency  factor  y.  Boundary 
layer  calculations  are  based  on  modified  equations  of  Patanker  and 
Spalding. 


STABLE  LAMINAR 
BOUNDARY  LAYER 
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Forest  uses  a  revised  pressure  gradient  parameter  m  =  mQV/v  =  f(Re^, 
y,  and  Tu) ,  which  accounts  for  the  eddy  viscosity  effect. 

Forest  found  fair  agreement  with  turbine  data  and,  like  Walker, 
found  large  areas  of  laminar  and  transitional  flow.  It  is  interesting 
tc  note  that  at  a  turbulence  level  of  6%,  transition  occurred  from  30% 
to  80%  over  the  chord  length  for  turbine  blades.  Forest's  model,  how¬ 
ever,  does  not  account  for  heat  transfer  effects,  which  Schlichting  con¬ 
siders  an  important  component.  Forest's  contribution  is  significant, 
since  he  accounts  for  the  Gortler  number  on  the  concave  surface. 

4.3.5.  Walker's  Correlation 

Walker  (8,  23)  presents  the  most  comprehensive  evaluation  of  tran¬ 
sition  on  compressor  blades.  Walker  conducted  tests  on  the  British  C4 

section  stator  blades  of  a  single  stage  research  compressor.  Chord 

4  5 

Reynolds  number  varied  from  2.5  x  10  to  1.9  x  10  ,  angle  of  incidence 

varied  from  -11°  to  +5°,  and  turbulence  intensity  ranged  from  2%  to  6%. 

Walker  observed  that  the  flow  was  basically  two-dimensional  over 

80%  of  the  stator  span  because  of  an  absence  of  secondary  or  radial  flows. 

He  then  used  a  stethoscopic  technique  and  baked-clay  tests  to  locate 

transition  experimentally.  He  analytically  located  the  instability 

point  Reg  crit  as  a  function  of  Stuart's  pressure  gradient  parameter 
2 

mg  *  6  /v(dU/dx).  He  used  a  mean  value  since  the  instability  point 
fluctuates  in  unsteady  flow.  Analytically,  transition  was  defined  as 
the  location  where  turbulent  spots  appear. 

Walker  also  observed  that  at  highly  negative  incidence  (approx  -6°) 
the  instability  region  occupied  45%  of  the  chord  length.  At  positive 
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angles  of  incidence  (approx  +4°),  the  instability  region  was  reduced  to 
20%  of  chord. 

Walker  detected  Tollmien-Schlichting  waves  in  the  trough  of  the 
instantaneous  "nsteady  velocity  trace,  and  observed  that  the  breakdown 
process  occurred  during  the  acceleration  phase.  At  the  unsteady  Reynolds 
number  of  ReNg  =  1000,  he  observed  that  transition  occurred  in  the 
"aperiodic"  mode  as  suggested  by  Obremski  and  Fejer  (28).  He  observed 
that  the  location  of  transition  increased  as  the  angle  of  incidence 
decreased,  but  that  transition  location  was  independent  of  Rec.  He 
also  noted  that,  in  general,  (Re^  t  -  Re^  decreased  as  Rec  decreased. 

Walker  then  correlated  the  location  of  transition,  based  on  the 
location  of  instability,  to  the  mean  shape  factor  Hffl  across  the  insta¬ 
bility  region.  Using  a  least  square  fit  to  Granville's  isolated  airfoil 
data  and  his  own.  Walker  obtained  the  following  transition  Reynolds 
number  correlation: 


(Refl  -  Re  .)/Refl  -  1.70  -  0.32  x  H 

o,tr  o,mean  m 


(28) 


4  7 

This  correlation  applies  to  3  x  10  <  Rec  <  5  x  10  .  Walker  also  pro¬ 

posed  a  transition  location  based  on  the  mean  shape  factor.  In  all  cases, 
the  mean  shape  factor  is  a  function  of  pressure  gradient  dp/dx  and  Re£ 
which  are  the  dominant  factors  affecting  amplification  rates. 

Walker's  correlation  does  not  include  freestream  turbulence,  for 
reasons  discussed  in  the  previous  section  on  freestream  turbulence 
effects.  Stated  briefly,  he  feels  that  turbulence  indirectly  affects 
transition  through  a  direct  modification  of  the  incidence  angle  and 


however,  were  not  varied  independently  of  Reynolds  number.  The  turbu¬ 
lence  intensity  decreased  with  decreasing  Reynolds  number.  This  may 
represent  a  flaw  in  Walker's  correlation. 

Finally,  Walker  observed  that  the  transition  length  x^,  -  x  is  nearly 
constant,  occu^  g  about  10%  to  20%  of  chord  for  negative  angles  of 
incidence.  Walker's  model  also  predicts  an  overly  large  xt  since  the 
shape  factor  H  is  too  large  in  his  calculation  scheme  (Thwaites'  method). 
Consequently,  laminar  separation  is  predicted  too  often  in  Walker's 
model. 

4.3.6.  The  Albers  and  Gregg  Method 

Albers  and  Gregg  (41)  modified  an  existing  finite  difference  boundary 
layer  calculation  computer  program  to  provide  an  improved  transition 
prediction  method.  The  instability  point  is  determined  from  Re^  crit  = 
f(mg,  St).  For  transitional  boundary  layer  computation,  a  turbulence 
model  is  employed  including  the  calculation  of  an  effective  viscosity 
and  intermittency  factor.  Transition  is  defined  as  the  location  where 
Y  =  0.01. 

To  determine  the  location  of  transition,  ReQ  and  ReQ  are  calcu¬ 
li  t)  ,tr 

lated  at  each  streamwise  station.  Calculation  of  ReQ  .  includes  the 

9, crit 

effects  of  pressure  gradient,  freestream  turbulence,  surface  roughness, 
suction,  and  surface  curvature.  Each  effect  is  normalized  into  a 
reduction  factor,  based  on  the  transition  Reynolds  numbers  for  each 
effect  presented  in  graphic  form  by  Schlichting  (5,  1960  edition). 

The  corrected  transition  Reynolds  number  is  computed  from  the  following 
equation 
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ReQ  =  Re  x  R,  x  R  Qx  R,,x  R„  x  R 
0,tr  ok  mQ  H  Tu  curv 

where  ReQ  represents  the  flat  plate  transition  Reynolds  number  and  R^  ^ 
represents  the  appropriate  reduction  factors.  This  calculation  is  done 
at  each  station. 

4.3.7.  The  McDonald  and  Fish  Method 

McDonald  and  Fish  present  a  finite  difference  computation  scheme 
for  calculating  transitional  boundary  layer  characteristics.  Transition 
is  located  at  the  point  of  minimum  skin  friction  coefficient. 

McDonald  and  Fish  employ  a  turbulence  model  approach  using  a 
turbulent  energy  equation  consisting  of  a  convection  integral,  a  net 
production  integral,  and  a  normal-stress  production  integral.  In  the 
initial  stage  of  transition,  the  normal-stress  production  term  equals 
zero.  The  velocity  distribution  and  boundary  layer  thicknesses  are 
then  computed  at  each  station.  McDonald  and  Fish  note  that  in  the  tran¬ 
sition  region,  the  viscous  stresses  and  apparent  turbulent  stresses  are 
approximately  equal. 

Freestream  turbulence  is  accounted  for  as  a  source  term  in  the 


turbulent  energy  equation.  This  changes  the  magnitude  of  the  Reynolds 
stress,  thus  decreasing  the  transition  region  length.  If  the  roughness 
is  increased  sufficiently,  a  roughness-induced  transition  begins. 
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5.  SUMMARY 

The  foregoing  review  should  indicate  the  complexity  of  transition 
prediction,  even  in  two-dimensional,  incompressible  flows.  The  inter¬ 
actions  of  the  parameters  are  not  fully  understood  or  universally  agreed 
upon;  the  interaction  of  freestream  turbulence  and  pressure  gradient  is 
only  one  example  of  this  disagreement.  Consequently,  data  correlations 
of  airfoil,  cascade,  and  compressor  boundary  layer  measurements  offer 
the  best  solution  to  transition  prediction.  Finite  difference  computa¬ 
tions  provide  a  means  for  calculating  transitional  boundary  layer  growth, 
but  still  rely  on  the  data  correlations  to  pinpoint  the  initiation  of 
transition. 

Stability  theory  provides  a  useful  guide  for  determining  the  overall 
effects  of  various  parameters,  even  though  it  accounts  for  only  one 
step  in  the  entire  transition  process.  Parametric  effects  other  than 
pressure  gradient  have  not  been  successfully  modeled. 

The  mechanism  of  flat  plate,  zero  pressure  gradient  transition 
appears  to  be  universally  accepted  thanks  to  the  large  number  of  flow 
visualization  tests  conducted.  However,  some  disagreement  exists  regard¬ 
ing  the  parametric  influences,  and  the  mechanism  of  turbulence  is  not 
agreed  upon  because  of  the  complexity  of  describing  its  characteristics. 

Airfoils  and  compressor  blades,  therefore,  present  a  difficult 
problem  since  all  of  the  parameters  mentioned  in  this  paper  act  simul¬ 
taneously.  Data  correlations  from  different  sources  usually  show  varia¬ 
tions  on  the  order  of  10  to  30%  in  the  transitional  Reynolds  number 
because  of  variety  in  blade  section  shapes  and  free  stream  disturbance 
spectra,  blade  roughness,  and  hidden  (or  unreported)  three-dimensional 
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effects.  This  author  suggests  that  the  second  derivative  of  the  pres- 
2  2 

sure  gradient  d  p/dx  may  help  account  for  the  blade  section  shape. 

Walker  (8),  however,  considers  the  shape  factor  H  capable  of  accounting 
for  this  effect. 

In  the  absence  of  a  fully  developed  transition  theory  (a  develop¬ 
ment  Morkovin  (2)  considers  improbable) ,  this  author  recommends  the  use 
of  Walker's  data  correlation.  Walker  has  presented  the  most  detailed 
study  of  the  transition  process  on  compressor  blades  and  has  also  sup¬ 
ported  many  of  his  findings  with  accepted  boundary  layer  theories,  par¬ 
ticularly  those  of  Schlichting  (5)  and  Obremski  and  Fejer  (28).  However, 
some  doubt  about  broad  application  of  Walker's  findings  remains  since 
he  did  not  vary  turbulence  intensity  independently  of  the  chord  Reynolds 
number.  However,  his  observations  of  freestream  disturbances  effects 
are  well  taken,  and  his  observations  regarding  laminar  separation  are 
in  agreement  with  the  observations  of  Roberts  (20)  and  Kiock  (22). 

Finally,  work  on  transition  is  essentially  incomplete  and,  as 
White  (1)  remarks,  "the  final  report  on  transition  may  never  be  handed 
in."  The  problem  remains  unresolved.  Consequently,  the  transition 
prediction  methods  summarized  in  this  paper  should  be  applied  cautiously 
and  with  an  awareness  of  the  fluid  mechanics  involved  in  each  particular 
flow  problem. 
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